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Application of a significantly large bias voltage to small Bi 2 Sr2CaCu 2 08 + , mesa structures leads to persis- 
tent doping of the mesas. Here we employ this effect for analysis of the doping dependence of the electronic 
spectra of Bi-2212 single crystals by means of intrinsic tunneling spectroscopy. We are able to controllably 
and reversibly change the doping state of the same single crystal from underdoped to overdoped state, without 
changing its chemical composition. It is observed that such physical doping is affecting superconductivity in 
Bi-2212 similar to chemical doping by oxygen impurities: with overdoping the critical temperature and the su- 
perconducting gap decrease, with underdoping the c-axis critical current rapidly decreases due to progressively 
more incoherent interlayer tunneling and the pseudogap rapidly increases, indicative for the presence of the 
critical doping point. We distinguish two main mechanisms of persistent electric doping: (i) even in voltage 
contribution, attributed to a charge transfer effect, and (ii) odd in voltage contribution, attributed to reordering 
of oxygen impurities. 

PACS numbers: 74.72.Gh, 74.62.Dh, 74.55.+V, 74.72.Kf 



I. INTRODUCTION 

High temperature superconductivity (HTSC) in cuprates 
occurs as a result of doping of a parent antiferromagnetic Mott 
insulator and properties of cuprates change significantly with 
doping. Superconductivity in overdoped cuprates is fairly 
well described by the conventional BCS-type second-order 
phase transition JTJO. But properties of underdoped cuprates 
are abnormal due to the persistence of the normal state pseu- 
dogap, strong superconducting fluctuations, or possibly pre- 
formed pairing, [4] and magnetism [5| at T > T c , There are 
indications that the transition from the normal to the abnormal 
behavior occurs abruptly at a critical doping point HHSHS). 
This may be a consequence of the quantum phase transition - 
a phase transition, which occurs at T — 0, in frustrated sys- 
tems as a result of a competition of coexisting order parame- 
ters. The coexistence of superconductivity at T < T c with the 
pseudogap lfT0lfl"2l . charge and spin density order Il8l l9l [T3l 
was indeed reported by several techniques. Clearly, detailed 
doping-dependent studies are needed both for understanding 
the puzzling nature of HTSC in cuprates and for the develop- 
ment of novel HTSC materials. 

Usually, the mobile carrier concentration is controlled by 
chemical doping via chemical substitution, or in case of 
cuprates also by variation of the oxygen content via appropri- 
ate annealing and subsequent quenching to room temperature. 
This allows an accurate control of the chemical composition, 
but less so of the local arrangement of impurities and disor- 
der, which is equally important for cuprates [ 14 1. For exam- 
ple, it is well established that properties of the YBa2Cu306 +x 
compound strongly depend not only on the average oxygen 
concentration but also crucially on the order/disorder of oxy- 
gen atoms in Cu-O-Cu chains fl31[T6ll . Therefore, analysis of 
the doping phase diagram of cuprates requires accurate con- 
trol of both, the concentration and the microscopic structure 
of impurities. 
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The carrier concentration can be also varied via two phys- 
ical doping processes, well established for semiconductors: 
photo-doping lfT7H T9l and through the electric-field effect 
E0l - l25l . In case of cuprates, physical doping may be per- 
sistent at low temperatures in a sense that it is relaxing very 
slowly after removing the light G7HHI or field IZTll25ti25l. 
Recently, a persistent electric doping via strong current in- 
jection was discovered 11241 [251 . It is resembling a resistive 
switching phenomenon in memristor devices |26 27] and is 
related to previous similar observations in point contact ex- 
periments on Bi2Sr2CaCu20g +A (Bi-2212) [23 1. Such an elec- 
tric doping is reversible, reproducible and easily controllable. 
It opens a possibility to analyze the doping dependence of 
HTSC on one sample without changing its chemical compo- 
sition [19]. Despite that, there were very few direct spectro- 
scopic studies of cuprates employing physical doping tech- 
niques. 

Intrinsic Tunneling Spectroscopy (ITS) provides an unique 
opportunity to probe bulk electronic properties of HTSC El . 
This technique utilizes weak interlayer (c-axis) coupling in 
quasi two dimensional HTSC compounds, in which mobile 
charge carriers are confined in Cu02 planes separated by some 
blocking layer (e.g. SrO-2BiO-SrO in case of Bi-2212). This 
leads to a formation of atomic scale intrinsic tunnel junctions, 
and to an appearance of the intrinsic Josephson effect at T < 
T c HUMS]. 

In this work we employ the persistent electric doping for 
analysis of the doping dependence of electronic spectra of 
Bi-2212 single crystals by means of intrinsic tunneling spec- 
troscopy l2l[3ll6l [T0l[32l . Controllable and reversible persis- 
tent physical doping is achieved by applying a c-axis voltage 
of a few volts to small Bi-2212 mesa structures. Thus we 
are able to change the doping state of Bi-2212 single crystals 
without changing its chemical composition. A wide doping 
range from a moderately underdoped to strongly overdoped 
state could be reached. It is observed that the physical doping 
is affecting the intrinsic tunneling characteristics of Bi-2212 
similar to chemical doping [6|. With overdoping the critical 
temperature and the superconducting gap decrease. With un- 
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derdoping the pseudogap rapidly increases, indicative for the 
presence of the quantum critical doping point in the phase dia- 
gram, and the c-axis critical current density rapidly decreases, 
indicating a progressively more incoherent interlayer tunnel- 
ing. We distinguish two main mechanisms of persistent elec- 
tric doping: (i) an even in voltage contribution, attributed to a 
charge transfer effect, and (ii) an odd in voltage contribution, 
attributed to reordering of oxygen impurities. 

The paper is organized as follows. In Sec. II we make a 
brief overview of physical doping mechanisms of cuprates. 
Sec. Ill provides experimental details. In Sec. IV we present 
the main experimental results and in Sec. V we discuss pos- 
sible mechanisms of persistent electric doping, followed by 
conclusions. 



II. PHYSICAL DOPING OF CUPRATES 

A. Photo-doping 

Photo-doping allows a wide-range variation of doping in 
the same sample lfT7lfl9l . The ordinary non-equilibrium 
photo-doping is quickly relaxing because of a very short life 
time (~ps) of photoinduced charge carriers Ifl9l . However, in 
underdoped YBa2Cu306+ x and some other cuprates a differ- 
ent type of persistent photo-doping takes place lfT7l [T8ll . It in- 
volves significant energies ~ 1 eV, which makes it metastable 
at low temperatures. Several mechanisms are contributing 
to the persistent photo-doping l33l . such as charge transfer, 
which changes the redox state of the impurity atom 1(171 . 
and ordering of oxygen impurities in the lattice Ifl8l . Photo- 
doping always leads to an increase of the doping level with 
respect to the initial state. 

B. Electric field effects 

Electric fields may both increase or decrease the number of 
mobile charge carriers, depending on the direction of the ap- 
plied field ll20l . The ordinary electric field effect is not persis- 
tent and exists only during the time an electric field is applied. 
Since the electric field penetrates only to the Thomas-Fermi 
charge screening length, A-^p < 1 nm, just a thin surface layer 
can be modified |34|. 

A persistent electrostatic field-effect due to net electric po- 
larization or trapped charges can be realized at the interface 
between a superconductor and a ferroelectric [ 35 1 or polar in- 
sulator [ 36 1 . This is also a surface phenomenon, but in case 
of layered cuprates, which represent stacks of metallic CuO 
planes sandwiched between polar-insulating layers [37 1, elec- 
trostatic charging of insulating layers may in principle lead to 
the bulk persistent electrostatic field-effect. 

Another type of a persistent and bulk electric -field effect 
has been observed at large current densities [21]. Similar to 
photo-doping, it was attributed to a charge transfer 11381 and 
reordering of oxygen impurities 11221 . Significant oxygen mo- 
bility in intense electric fields also leads to a resistive switch- 
ing phenomenon ||23| . 



C. Resistive switching in complex oxides 

The resistive switching phenomenon occurs in many com- 
plex oxides and is the basis for the development of resistive 
memory devices. Several mechanisms may be involved in the 
resistive switching phenomenon 112611271 . including a change 
of the redox state of some of the elements, oxygen migra- 
tion, and filament formation. Resistive switching has been 
observed on depleted surfaces of Bi-2212 cuprates 11231 and 
attributed to oxygen migration. Recently, it was demonstrated 
that the resistive switching-like phenomenon can be used for 
controllable and reversible doping of small Bi-2212 micro- 
structures over a wide doping range l24ll25l . 



HI. EXPERIMENTAL 

Mesas were fabricated on top of freshly cleaved Bi-2212 
single crystals by means of optical lithography, Ar ion milling 
and focused ion beam trimming. Four batches of crys- 
tals were used: pure near optimally doped (OP) Bi-2212, 
pure strongly underdoped (UD) Bi-2212, lead-substituted 
Bii.75Pbo.25Sr 2 CaCu20 8+I 5 [Bi(Pb)-2212], and yttrium-doped 
Bi 2 Sr 2 Cai_ x Y x Cu 2 8+5 [Bi(Y)-2212]. Mesas of different 
sizes from 5 X 5 to 1 x 0.5 ^im 2 and with a different number of 
junctions N = 8-56 were studied. Details of the sample fabri- 
cation can be found in Ref. [ 39 1 . All studied mesas exhibited a 
persistent doping effect upon application of a sufficiently large 
bias voltage. 

The samples were placed in a flowing gas cryostat and mea- 
sured in a three-probe configuration with a common top gold 
contact. The ground contacts for current and voltage were 
provided through other mesas on the same crystal. A Keithley 
K6221 current source and a FPGA-based arbitrary waveform 
generator and lock-in amplifier were used to bias and measure 
the samples. Biasing was done at pseudo-constant voltage, 
with an optional small superimposed ac voltage to simultane- 
ously measure the high-bias differential resistance in addition 
to the dc resistance. Positive bias is defined as current going 
into the mesa through the common top contact, as sketched in 
Fig. [T](a). 

IV. RESULTS 

From chemical (oxygen) doping studies it is known that 
doping/undoping of Bi-2212 is accompanied by a propor- 
tional decrease/increase of the c-axis resistivity [6|. There- 
fore, we can control the doping state by tracing the mesa re- 
sistance. 



A. Dynamics of electric doping 

The basic features of the dynamics of the persistent electric 
doping are shown in Fig.[T] Panel (a) shows the time-evolution 
of the Bi(Pb)-2212 mesa resistance at a bias of »1.9V at 
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FIG. 1. (Color online). Dynamics of electric doping (a-c) for a Bi(Pb)-2212 mesa at T = 135 K and (d-f) for an optimally doped Bi-2212 mesa 
at T = 2K. Panel (a) shows time evolution of the dc resistance at a bias of ss 1.9 V. The dashed line represents a stretched exponential decay. 
The inset schematically shows a sketch of a mesa structure and the direction of electric field at positive bias, (b) Demonstration of odd-in- 
voltage doping: the mesa resistance increases/decreases at positive/negative bias. Panel (c) shows the predominantly odd voltage dependence 
of the logarithmic rate of the dc resistance change d In R ic /dt. Different symbols represent different runs. The inset in (c) demonstrates the sign 
change of the doping direction at higher bias. Panel (d) demonstrates a gradual change from a negative to positive doping rate with increasing 
bias voltage and time. It also demonstrates a history dependence of the doping rate, i.e., a different sign of resistance change at the same bias 
voltage, depending on the former doping treatment. The history dependence upon several sweeps of the bias voltage are shown in (e). Panel 
(f) demonstrates the predominantly even in voltage doping for OP Bi-2212. However an asymmetry indicates the presence of a subdominant 
odd in voltage doping. Note that despite a significant difference between the two mesas, the threshold doping voltage is similar ~ ±1.7 V (c). 



T — 135 K. It is seen that the mesa resistance is decreas- 
ing with time, indicating a gradual doping of the mesa. The 
doping rate decreases with time, following a stretched expo- 
nential decay, R — Rq + exp[-(t/T~f], shown by the dashed 
line in in Fig. [T|(a), which is also typical for persistent photo- 
doping 1 17, 18 1 . The doping can be equally well performed 
at any temperature from 4 to 300 K, but the rate is increasing 
with T. In most cases we perform doping at low T in order to 
be able to immediately probe the superconducting character- 
istics. Upon reduction of the bias below the threshold voltage 
the state of the mesa remains stable even at room temperature 
on the time scale of several days. 

The resistive change is reversed upon voltage reversal, as il- 
lustrated in Figs.[T](b) and (c) for the same Bi(Pb)-2212 mesa. 
The resistance decreases for positive bias (electric field into 
the crystal) and increases for negative bias (electric field to- 
wards the top contact), which indicates that we can control- 
lably and reversibly dope and undope the mesa. The doping 



rate and direction depend both on the sign and the absolute 
value of bias voltage. 

Figure [T](c) summarizes the bias dependence of the doping 
rate for the Bi(Pb)-2212 mesa. Below the threshold voltages, 
l^dcl ^ 1.7 V, the mesa resistance is stable. Upon increasing 
the bias voltage, the resistance of the mesa starts to gradually 
change at a rate that increases drastically up to |V| ~ 2.2 V 
as shown in the inset of Fig.[T](c). A further voltage increase 
reduces the rate and then reverses the resistance alteration rate 
(see the inset at V = -2.32 V). The sign change of the alter- 
ation rate at high bias is in agreement with the observations by 
Koval etal. [24|. The behavior in this regime is, however, his- 
tory dependent, as may be seen from Fig. [T](c), and the final 
state depends on how long time the mesas was biased at every 
bias voltage. The doping process for the Bi(Pb)-2212 mesa, 
Fig.[T](c), is predominantly odd in bias voltage, i.e., the direc- 
tion of doping is changed when the sign of the bias voltage is 
changed. 
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FIG. 2. (Color online) Intrinsic tunneling characteristics of the same OP Bi-2212 (a-d) and Bi(Pb)-2212 (e-g) mesas as in Fig. 1 at different 
doping states, (a) T dependence of the ac resistance R a{: (T) in the initial high-resistive state (1) and subsequent low resistance doping states 
(4-6). The inset illustrates the significant enhancement of T c of the depleted surface junction. Panels (b-d) represent different parts of / - V 
characteristics in the doping states 1, 2, 3 and 6 at T = 2 K. A variation of the tunnel resistance, sum-gap kink and the critical current is clearly 
seen. Panels (e) and (f) show R X (T) for Bi(Pb)-2212 in the initial high-resistive state (1) and subsequent low resistance doping states (4) and 
(6). This mesa with a large number of junctions N = 56 exhibits a significant spread in the T c of individual junctions, seen as small resistance 
drops in (f). Panel (g) shows / - Vs of the Bi(Pb)-2212 mesa in the doping states 2, 4 and 5. A progressive increase of the critical current and 
decrease of the sum-gap kink voltage is seen. 



Figures [T](d) and (e) show a detailed view of the time and 
voltage dependence of doping for a near optimally doped pure 
Bi-2212 mesa. The top panel of (d) shows the time evolu- 
tion of the resistance of an OP Bi-2212 mesa for different bias 
voltages, shown in the bottom panel. It is seen that the resis- 
tance is constant at V — 1 .7 V, and starts to increase slowly 
at 1.9 V. However, at 2 V the resistance initially increases but 
then starts to decrease after a few minutes. This clearly shows 
that there are two counteracting processes: a positive rate 
mechanism that saturates quickly and a mechanism with nega- 
tive rates that is dominating at longer times and at higher volt- 
ages. The second process also saturates with time, which is 
clear from Fig. [T](a) and history dependent rates of Fig. [T](e). 
At larger voltage, the resistance steadily decreases at a rate 
which is strongly bias dependent, as shown in Fig.[T](f). 

Figures [T](e,f) show the bias dependence of doping and a 
history dependence upon sequential voltage sweeps (e). It is 
seen that for the OP Bi-2212 mesa the electric doping is pre- 
dominantly even in voltage, i.e., the direction of doping does 
not depend on the sign of the bias voltage. However, cer- 
tain asymmetry of the doping rate versus bias voltage charac- 



teristics in Figs. 2 (c) and (f) indicates the presence of sub- 
dominant even- and odd-in-voltage contributions for Bi(Pb)- 
2212 and OP Bi-2212 mesas, respectively. 



B. Doping dependence of ITS characteristics 

Using the described method, the superconducting proper- 
ties of the mesas have been altered to different intermediate 
doping states denoted by a successive number. The electric 
doping changes all mesa characteristics: the c-axis resistivity, 
the critical temperature T c , the c-axis critical current density 
J c , the superconducting energy gap, A, the c-axis pseudogap 
and the c-axis resistivity in a manner very similar to chemical 
(oxygen) doping J6). 

Figure [2] shows temperature dependencies of low bias ac- 
resistances for the initial, high resistance state (HRS), and 
doped, lower resistance states (LRS), for (a) an OP Bi-2212 
and (e,f) a Bi(Pb)-2212 mesa. The I-V characteristics of those 
mesas at different doping states are presented in panels (b-d) 
and (g), respectively. 



5 




l/(V) 1/(V) 7"(K) 



FIG. 3. (Color online). Doping of a strongly underdoped Bi-2212 mesa, (a) I-V characteristics in the initial state 1 (black) and two successive 
doping states 2 (red) and 3 (green). The increase of the critical current by a factor five is seen, (b) The corresponding d/-dV characteristics. 
The superconducting sum-gap peak moves to slightly lower voltages and becomes sharper with doping. The c-axis pseudogap hump voltage 
rapidly decreases with doping, (c) Zero bias ac resistivity R^T). The T c has increased by about 15 K in the state 3. 



From Fig. |2](a) it is seen that the initial state was charac- 
terized by the strong thermal-activation-type increase of re- 
sistance with decreasing T, typical for underdoped Bi-2212 
QUI . The general shape of the resistive transition was de- 
scribed in Ref. J2) . At T c ~ 82 K the resistance dropped to the 
top-contact resistance, which originates from the first deteri- 
orated junction between the top CuO plane, shortly exposed 
to atmosphere after cleavage, and the second, un-deteriorated 
CuO plane [2|. Initially this junction had a very low T' c and 
a very small critical current, 7 C , as can be seen from the cor- 
responding I-V in panel (d). After electric doping, the resis- 
tance in the normal state dropped almost three times and be- 
came less semiconducting. The main T c of the mesa changed 
only slightly, indicating that the doping was changing around 
the optimal doping level with the fiat T c vs. doping depen- 
dence. However, the properties of the top junction changed 
drastically: the T' c increased to ~ 50 K, and the critical current 
increased ten-fold as shown in panel (d), even though it still 
remains ~ 20 times smaller than for the rest of the junctions 
in the mesa, as can be seen from panel (c). This indicates 
that the surface CuO plane was initially strongly underdoped 
and, therefore, responded much stronger to variation of dop- 
ing, due to the steep T c vs. doping dependence at the under- 
doped side of the doping phase diagram of cuprates J6). A 
similar trend was also observed in photo-doping |[T7l[T8l . 

Another possible reason for the stronger response of the un- 
derdoped top junction is the larger c-axis resistivity of under- 
doped intrinsic Josephson junctions [6|. Because of that the 
electric field is not uniformly distributed along the mesa but 
is larger in the high-resistive top junction. This together with 
the strong voltage dependence of the electric doping leads to a 
faster doping of the top junction and the doping may even go 
in the opposite direction with respect to the rest of the mesa. 

The effect of non-uniform doping along the height of the 
mesa becomes more pronounced in higher mesas with a larger 
number of intrinsic Josephson junctions. This is seen from 
R(T) for the Bi(Pb)-2212 mesa from Fig. |2](e), which con- 
tained a fairly large number of junctions N ~ 56. It is seen that 



after doping some junctions retained the initial T c ~ 90 K, but 
some were very strongly overdoped to T c ~ 30 K. Fig.[2](f) 
shows R(T) at the intermediate doping state 4. Small drops 
represent critical temperatures of individual junctions in the 
mesa. Apparently there is a gradual distribution of T c along 
the height of the mesa. 

Doping of Bi-2212 leads to a rapid increase of the c-axis 
critical current density (6). This is clearly seen from I-V 
curves for the Bi-2212 (OP) mesa shown in Figs.|2](b-d). The 
increase of 7 C of the surface junction is shown in panel (d). 
One-by-one switching of the rest of the junctions into the re- 
sistive state at I > I c leads to the appearance of multiple- 
quasiparticle (QP) branch structures in the I-V curves. The 
corresponding critical current at the first QP branch is shown 
in panel (c). It also strongly increases with doping. The same 
effect doubles I c in the Bi(Pb)-2212 mesa, in Fig. |2](g). 

From Fig.|2](b) it is seen that the I-V curves exhibit a kink at 
large bias, followed by an ohmic tunnel resistance. The kink 
represents the sum-gap singularity in superconducting tunnel 
junctions at V = 2A/e per junction [2 3|. This is the basis of 
the ITS technique, which allows analysis of the superconduct- 
ing energy gap A in the bulk of the Bi-2212 single crystal. 

Accurate analysis of the electronic spectra with the ITS 
technique requires mesas with a small area and a small num- 
ber of identical junctions. This is needed for avoiding pos- 
sible artifacts, associated with self-heating, in-plane non- 
equipotentiality, and spread in junction parameters l40l . 
This is particularly important for the analysis of the genuine- 
shape of tunneling characteristics, which remains a contro- 
versial issue PP . Even though the I-V characteristics in 
Figs. |2](b) and (g) are distorted by self-heating, as evident 
from a back-bending at large bias, the general trend for a 
variation of the sum-gap kink with doping is clearly seen: 
The superconducting gap decreases and the sum-gap kink be- 
comes sharper with (over-)doping. This qualitative conclusion 
is not affected by self-heating because the dissipation power 
at the kink decreases with decreasing resistance and becomes 
smaller with subsequent doping. Thus, with over-doping the 
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FIG. 4. (Color online). Intrinsic tunneling characteristics of a small Bi(Y)-2212 mesas at low T before (slightly underdoped, HRS) and after 
(optimally doped, LRS) doping by a short voltage pulse. Panels (a) and (b) show I-V curves at a large scale and at the quasiparticle branches, 
respectively. A reduction of the high bias tunnel resistance and a simultaneous increase of the critical current in the doped state is seen. Panel 
(c) shows dl/dV ITS characteristics. A change in the shape of the curves is seen: the peak becomes sharper and the dip-hump less pronounced 
in the optimally doped state, compared to the initial underdoped state. 



superconducting sum-gap singularity becomes sharper and 
moves to lower voltages despite the progressive reduction of 
self-heating. This clearly reveals the doping -variation of the 
genuine c-axis tunneling characteristics [41 1. A similar ten- 
dency was observed by other techniques, including the angu- 
lar resolved photoemission spectroscopy Il42l . scanning tun- 
neling spectroscopy [43 1, and tunneling spectroscopy on point 
contacts 11441 . as well as in previous ITS studies involving 
chemical (oxygen) doping |2] |6l [32] . 

Figure [3] shows the electric doping of an initially strongly 
underdoped Bi-2212 mesa. A five-fold increase in critical 
current in the I-V characteristics (a), and an increase of T c 
of about 15 K (c) can be seen. The tunneling conductance 
d//dV(V) curves are shown in panel (b). It is seen that the 
sum-gap peak shifts to slightly lower voltages and becomes 
sharper with doping. The hump voltage, attributed to the c- 
axis pseudogap [2|, rapidly decreases with increasing conduc- 
tance, pointing towards an abrupt opening of the pseudogap 
at the critical doping point, consistent with chemical doping 

studies 013 na. 

Note that in all studied cases the electric doping has lead to 
a significant increase of the critical current, while the super- 
conducting gap was decreasing. The anticorrelation between 
I c R n and A in underdoped Bi-2212 has been reported before 
[6 1 and was attributed to progressively more incoherent c-axis 
transport in combination with the li-wave symmetry of the su- 
perconducting order parameter. 



C. Short-pulse doping 

So far we were discussing a gradual electric doping of Bi- 
2212 mesas at the time scale of an hour, as shown in Fig. [T] 
Such a long time doping allows a very strong variation of the 
doping state, but often leads to inhomogeneous doping within 
the mesa height, as shown in Fig. [2] Koval et.al. 11251 demon- 
strated that a short-pulse doping strategy leads to highly re- 



versible and reproducible doping, similar to resistive switch- 
ing in point contacts [23 1. This is probably related to the lack 
of significant electromigration during the short pulse, which 
may eventually lead to an irreversible destruction of the crys- 
tal structure ED . 

Figure [4]represents the ITS characteristics at T ~ 30 K for 
a small ~ 2 x 2|int 2 Bi(Y)-2212 mesa with a small amount 
of junctions and small self-heating [40 1. In the initial high- 
resistive state the mesa is slightly underdoped with T c ~ 91 K. 
The mesa was switched to a low-resistive state by a short volt- 
age pulse V > 2 V of about a milliseconds width. The pe- 
riodicity of QP branches in Fig. |4](b) demonstrates that af- 
ter switching into the LRS the mesa remains highly uniform. 
From comparison of I-V curves in panels (a) and (b) it is seen 
that the decrease of resistance by ~ 1/3 is accompanied by an 
almost four-fold increase of 7 C . The T c increased to ~ 93 K be- 
ing an indication that the mesa became near optimally doped. 

Fig.|4](c) represents the tunneling conductance d//dV in the 
high-resistive (underdoped) and the low-resistive (optimally 
doped) states. The following main changes in ITS spectra are 
seen: the superconducting sum-gap peak voltage decreased 
in the doped low resistive state and the shapes of spectra are 
changed. The relative sum-gap peak height dI/dV(V p )R n is 
increased by about 50% in the low-resistive state. The high- 
resistive state exhibits a peak-dip-hump structure, which is 
less obvious in the low-resistive state. All this is similar 
to the slow-doping case, Figs. [2] and 3, and consistent with 
the change of doping from the slightly underdoped to near 
optimally -doped state [6|, in accordance with other spectro- 
scopic studies involving chemical (oxygen) doping ||42T|441 . 

As already discussed above, the increase of sharpness of 
the sum-gap kink in the I-V curve and the peak in the d//dV 
curve in the low-resistive state, reported in Figs. [2] and [4] and 
the change of the shape of the peak-dip-hump feature in Fig. [4] 
(c) can not be attributed to self-heating, because the actual 
power dissipation at the peak is decreasing in the low -resistive 
state. For example, the corresponding powers at the peaks in 
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FIG. 5. (Color online) Short pulse resistive switching in a strongly underdoped Bi-2212 mesa, (a) I-V characteristics in the initial high-resistive 
state (HRS) and the doped low-resistive state (LRS) at T = 2 K. Panels (b) and (c) demonstrate the resistive switching sequence _R ac versus 
time at T = 100 K. The switching was made by a train of short positive and negative pulses, shown in the bottom of panel (b). It is seen that 
the HRS is stable, but the LRS is initially relaxing and then saturates at a resistance below the HRS. Several hundred reproducible resistive 
switching events can be achieved without visible degradation (c). 



Fig. |4](c) are P = 0.21 and 0.17 mW, respectively. This ob- 
servation supports the conclusion of Ref. [41 1 that the appear- 
ance and the shape of the peak-dip-hump structure in the ITS 
characteristics of small mesas is determined primarily by the 
doping level. 

A detailed analysis of short pulsed resistive switching has 
been performed on strongly underdoped Bi-2212 mesas. For 
negative pulses with 50 [is length and compliance voltages up 
to -2.5 V there was no resistive switching. At -3.0 V, a re- 
duction of the quasiparticle resistance to 99.3% is observed 
(LRS). Nine subsequent pulses reduce the resistance further 
to 98.4%. A single positive pulse switches the resistance back 
to the initial HRS, while more positive pulses do not result 
in additional changes. A single negative pulse with a higher 
compliance voltage of -3.5 V instead reduces the resistance 
of the LRS to about 93.4% of the HRS, and the corresponding 
positive pulses switches it back. Pulses with lower compli- 
ance voltage lead to a partial switching but additional pulses 
with the same compliance do not lead to a significant change. 
Fig. |5](a) shows I-V curves of another UD Bi-2212 mesa at 
T = 2 K in the HRS and the LRS obtained with a 3.5 V pulse 
of 50 [is width. The general difference between I-Vs after the 
pulsed doping in Figs. |4](a) and |5](a) is the same as for the 
slow doping in Figs.[2](a) and 3 (a). 

Fig. |5](b) demonstrates a reproducible resistive switching 
between HRS and LRS at elevated T = 100 K for another UD 
Bi-2212 mesa. The switching was performed using a simi- 
lar positive and negative pulse sequence with +3.5 V, a pulse 
width of 100 [is and an interval between pulses of 3 s, shown 
in the bottom panel of Fig. |5](b). Panels (b) and (c) show the 
corresponding time sequence of the measured zero-bias resis- 
tance of the mesa. It is seen that negative voltage pulses lead 
to switching into the LRS while a subsequent positive pulse 
switches the mesa back into the HRS. The corresponding re- 
sistance change rates are of the order of d In R/dt » +1 000 s~ 1 , 
which are somewhat higher than for the slow doping shown in 
Fig. [TJ(d), but not inconsistent with that data, taking into ac- 



count that the compliance voltage is also significantly higher. 
It is seen that the HRS is stable and shows no visible relax- 
ation at T = 100K, while the LRS is initially relaxing with 
the characteristic time tlrs = (0.77 ± 0.30) s and then sat- 
urates before reaching the HRS, as shown in panel (b). At 
T = 280 K, the behavior is similar with a time constant of 
tlrs = (0.75 + 0.34)s. 



V. DISCUSSION: MECHANISMS OF PERSISTENT 
ELECTRIC DOPING 

In Sec. II we briefly reviewed known mechanisms of persis- 
tent physical doping of cuprates. It is likely that some of them 
are playing a role in persistent electric doping, studied here. 
Indeed, the phenomenon is clearly related to the persistent 
electric -field effect 11211 . observed in YBa2Cu3C»6 + x, which in 
turn is clearly related to persistent photo-doping observed for 
various cuprates lfT71[T8l . 

To identify possible mechanisms, we first summarize char- 
acteristic features of the persistent electric doping: 

i) Observed different voltage dependencies (odd and even) 
indicate that several distinct mechanisms are involved. 

ii) The doping rate shows a threshold-like behavior as a 
function of bias voltage (see Fig. [T]>. The threshold voltage de- 
pends on T in a thermal-activation manner, i.e. decreases with 
increasing T. Remarkably, the threshold voltage is weakly 
dependent on the number of junctions in the mesa, consis- 
tent with previous reports l24ll25l . and is comparable to that 
for a single point contact (231. This suggests that the phe- 
nomenon is connected to some characteristic energy ~ 1 eV, 
rather than directly to the electric field. Indeed, for a given 
voltage, the latter should scale inversely proportional to the 
number of junctions in the mesa, i.e., would not be universal 
for different mesas (that is why we hesitate to refer to the phe- 
nomenon as a persistent electric field effect, and rather call it 
persistent electric doping). 
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iii) However, the role of the electric field should not be 
underestimated. Indeed, the displacement field is given by 
D = VsJNt, where s T is the dielectric constant and t is the 
thickness of the insulating barrier between CuO bi-layers. The 
ratio t/s r ^0.1 nm was estimated from an analysis of Fiske 
(geometrical resonance) step voltages OTI . Therefore, the dis- 
placement field in intrinsic junctions is 

D^-xl0 8 cm _1 . (1) 
N 

For V s 2 V and N - 10, which corresponds to the case of 
Fig. [4j one would get a very large value D 2 x 10 7 V/cm, 
which is certainly capable of seriously polarizing and displac- 
ing ions in complex oxides ll2Tll23ll26ll . 

iv) The phenomenon is not associated with a net change of 
the oxygen content, which may only decrease within the cryo- 
stat. To the contrary, mesas can be repeatedly and reversibly 
doped and undoped, as shown in Figs.[T](b) and |5](b). 

A. Charge transfer and electrostatic charge trapping 

If an injected electron has a high enough energy, it may join 
one of the ions, leading to a change of the redox state and the 
effective doping. It was suggested that such a "charge trans- 
fer" mechanism is involved both in photo-doping of cuprates 
IfTTl and in the resistive switching phenomenon in other com- 
plex oxides [26|. 

Alternatively, the electron may be trapped (localized) in di- 
electric parts, leading to electrostatic charging of the sample. 
The Bi-2212 compound has a layered structure with metallic 
CuO planes sandwiched between polar insulator BiO layers 
OTI . In this case the electrostatic charging will take place in 
BiO layers, which may affect the doping state of the neighbor- 
ing CuO planes via the electrostatic field-effect. The electro- 
static charging of means by current injection takes place uni- 
formly within the whole structure. Therefore, unlike the con- 
ventional electric field effect ll20ll and the electrostatic field ef- 
fect at the interface between a superconductor and a ferroelec- 
tric material Il35ll36l . the current injection may lead to a per- 
sistent bulk electrostatic field-effect doping of Bi-2212. Koval 
et. al. [24] emphasized the similarity of the phenomenon with 
the floating-gate effect utilized in Flash memory devices. 

Both types of charging effects have common similarities: 

i) The charge transfer requires a certain energy (~eV for 
YBa2Cu306+xX rather than electric field. 

ii) The sign of the current and the direction of the electric 
field does not matter. Therefore, such doping should be even 
with respect to the voltage sign, consistent with the observa- 
tions in Ref. Il24l . 

Therefore, we attribute even in voltage persistent electric 
doping to charge transfer and/or charge trapping mechanisms. 

B. Oxygen reorientation and reordering 

It is well established that the doping state of cuprates de- 
pends not only on the amount of off-stoichiometric oxygen, 



but also on the relative orientation of the oxygen bonds ifTSl . 
Therefore, the doping state can be changed by oxygen reorder- 
ing. Since the required energy is large ~eV, compared to ther- 
mal energies, oxygen reordering is a slow process and does 
not take place spontaneously at low enough T. Oxygen re- 
ordering is considered as one of the main mechanisms of the 
persistent photo [ 18 1 and electric field ETll doping. 

In the case of persistent electric-field doping, the oxygen 
reordering is steered by the polarization. Therefore, the di- 
rection of doping should depend on the direction of the elec- 
tric field, i.e., should be odd with respect to the bias volt- 
age. We clearly see such a contribution in our experiment, 
see Fig. [T|(b). Note that the odd-in-voltage contribution was 
reported in the point-contact case [23 1, but not reported in pre- 
vious related works made on zig-zag type Bi-2212 microstruc- 
tures [24 25]. The geometry of the latter samples is symmet- 
ric with respect to the electric field direction (changing the 
sign of the electric field is equivalent to flipping their sam- 
ple up-side down). This is not the case in point contacts and 
mesa structures, studied here, for which the fields down (into 
the crystal) and up (into the top electrode) are not equivalent. 
Therefore, the difference may partly be due to the difference 
in sample geometry, or to the observed sample-dependence of 
the relative strength of odd and even in voltage doping contri- 
butions, as shown in Figs. 1 (d) and (e). 

Thus, we attribute the odd in voltage persistent elec- 
tric doping mechanism to field-induced oxygen reorienta- 
tion/reordering. 



C. Irreversible processes: electromigration, filament and arc 
formation 

An increase of the bias voltage above ~ 2.5 - 3.0 V leads 
to a gradual increase of the current and irreversible change 
of the mesa properties. A similar phenomenon was observed 
in YBa2Cu306 +x thin films and attributed to electromigration 
and field-induced diffusion of oxygen, which is even in bias 
voltage. The increase of conductance is probably due to a 
dielectric breakthrough in the insulating BiO layers, which 
leads to a pin-hole and filament formation. Thus, we at- 
tribute the slow and irreversible drift of the mesa character- 
istics at large bias voltages to electromigration in the mesas. 
This destructive process is, however, distinctly different from 
the reversible and reproducible electric doping effect, reported 
above. 

After deterioration by electromigration, the mesa character- 
istics become similar to resistive switching characteristics for 
point contacts on top of oxygen-depleted, Bi-2212 surfaces 
ll23l . At even higher bias the resistance becomes very high 
(infinite). But an inspection in a microscope shows that the 
mesa itself remains intact. There is no physical evaporation of 
material or a crater at the place of the mesa, as in the case of 
a violent electric discharge. Instead there are clear indications 
of an arc formation at one of the sharp corners of the mesa, 
which probably leads to delamination of the structure and me- 
chanical disattachment of the mesa from the base crystal. 
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D. Mechanisms of energy accumulation 

The most puzzling property of the reported persistent elec- 
tric doping is that the required bias voltage is weakly depen- 
dent on the number of junctions in the mesa Il24l [25). This 
is clearly seen from the presented data, for which the thresh- 
old voltage is always ~ 2 V for N = 9 in Fig. [4] N = 1 1 in 
Fig.[T](f) and N ~ 56 in Fig.[T](d), which is also similar to that 
for a single point-contact ~ 1.5 V [23). It is, therefore, clear 
that electric doping requires a certain electron energy, rather 
than electric field. However, for stacked tunnel junctions, the 
energy acquired by the injected electron in every tunneling 
event is proportional to the voltage drop across the junction, 
5E st eV/N, and is significantly smaller than the threshold 
energy. The main question is, therefore, how the electrons 
accumulate a sufficiently large energy, required for doping. 

In Ref. 1125) it was suggested that an electron can accumu- 
late energy upon sequential tunneling through several junc- 
tions without relaxation. However, in this case only electrons 
in the last junction will have enough energy. This would result 
in a strongly inhomogeneous doping in different junctions. 
Moreover, the probability of sequential tunneling without re- 
laxation is small, because of the small ratio of relaxation time 
(t ~ ps) [45 1 to the tunnel time, f tun , t tun /r » 1. The probabil- 
ity of sequential tunneling through Adjunctions is decreasing 
rapidly oc (T/t tun ) N with increasing N. This should lead to a 
dramatic increase of the doping time with increasing N. In- 
deed, suppose that it takes fa ~ lOmin for Adjunctions at 
7=1 mA to dope the mesa. This will involve N e - Ifa/e tun- 
neling events in each junction (e is the electron charge). Since 
the probability of sequential tunneling through 2N junctions is 
decreasing quadratic ally, it would require tunneling events 
per junction, which will take tzN = N e t^ = 3.75 x 10 19 min. 
However, such a dramatic increase of the doping time with 
increasing mesa height is inconsistent with experiments. 

For the sequential tunneling scenario to be relevant, the ra- 
tio ?tun/i" should rapidly drop with increasing electron energy 
and become of the order of unity at E ~ 1 eV. This may be 
caused by resonant tunneling, which increases the tunneling 
rate of quasiparticles with a certain energy, and/or by a dras- 
tic slowing down of the high energy quasiparticle relaxation, 
which may be caused by a rapid decrease of the Eliashberg's 
electron-boson spectral function and a gap in the correspond- 
ing bosonic density of states at high energies B31l . Such a 
scenario is interesting to investigate because it may give an 
information about the bosonic spectrum, involved in Cooper 
pairing, and thus provide a clue about the electron-boson cou- 
pling mechanism, responsible for high-T^ superconductivity 
in cuprates [45 1. 

We also want to propose an alternative mechanism for the 
energy accumulation of electrons: the formation of electric- 
field domains in the natural atomic superlattice formed by the 
mesa. Electric field domains are well studied in semiconduct- 
ing superlattices |46|. They appear in weakly coupled super- 
lattices close to the resonant tunneling condition. The corre- 
sponding non-linearity leads to an instability and a multiple- 
valued current-voltage characteristics. As a result, the electric 
field distribution in the superlattice becomes nonuniform and 



is concentrated in one or several junctions. 

A possibility of a formation of electric field domains in Bi- 
2212 mesas is not just a hypothesis. In fact, the multiple- 
branch I-V of Bi-2212 mesas due to one-by-one switching of 
intrinsic junctions from the superconducting to the resistive 
state, shown in Fig.|4](b), is due to a formation of electric field 
domains in individual tunnel junctions. The formation of elec- 
tric field domains in Bi-2212 mesas at high bias would explain 
many of the features of the studied persistent electric doping. 
In this case electrons in the domain may gain an energy close 
to eV without sequential tunneling through the whole mesa. 
Furthermore, since domains are typically dynamic and prop- 
agate through the whole superlattice [46], this would also ex- 
plain the uniformity of doping in the whole mesa, and not just 
in the outermost junction. 



VI. CONCLUSIONS 

We have studied the effect of persistent electric doping 
on intrinsic tunneling characteristics of small Bi-2212 mesa 
structures. It was shown that the application of a sufficiently 
large voltage to the mesas leads to a controllable and re- 
versible physical doping of the mesas, without a modification 
of their chemical composition. This allows the analysis of 
bulk electronic spectra in Bi-2212 in a wide doping range on 
one and the same mesa. This physical doping has the same 
effect as chemical (oxygen) doping on the intrinsic tunneling 
characteristics of Bi-2212: the c-axis resistivity decreases, the 
critical current increases and the energy gap is decreasing to- 
gether with T c with over-doping. The anticorrelation between 
I c R n and A indicates that the c-axis transport becomes progres- 
sively more incoherent at the underdoped side of the phase 
diagram. An analysis of the doping variation of the intrin- 
sic tunneling characteristics of the same mesa provides a clue 
about its genuine shape: with subsequent doping, the sum- 
gap peak in the tunneling conductance becomes sharper and 
the pseudogap hump rapidly decreases with doping, suggest- 
ing the presence of a critical doping point, in agreement with 
previous chemical doping studies (6). 

By analyzing the bias- and time dependence we could iden- 
tify different mechanisms involved in the persistent electric 
doping: i) The even-in-voltage process via charge transfer 
and/or charge trapping, ii) The odd-in-voltage process via 
oxygen reordering. Those are distinct from the irreversible 
electromigration and oxygen electrodiffusion, observed at 
higher bias. 

We confirm the previous report ll25ll that the threshold volt- 
age for the electric doping is weakly dependent of the number 
of junctions in the mesas and is similar to that for a single 
surface point contact ll23l . This indicates that it is the en- 
ergy of injected electrons, rather than electric field, that de- 
termines the phenomenon. We suggest that the required en- 
ergy accumulation by tunnel electrons may be due to a forma- 
tion of electric field domains in the natural atomic superlattice 
formed in the Bi-2212 single crystal. 
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